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ABSTRACT 

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we have detected emission 
in the S(0), S(l), and S(2) pure-rotational (u = — 0) transitions of molecular hydrogen (H2) towards 
6 positions in two translucent high Galactic latitude clouds, DCld 300.2—16.9 and LDN 1780. The de- 
tection of these lines raises important questions regarding the physical conditions inside low-extinction 
clouds that are far from ultraviolet radiation sources. The ratio between the S(2) flux and the flux 
from PAHs at 7.9 microns averages 0.007 for these 6 positions. This is a factor of about 4 higher than 
the same ratio measured towards the central regions of non-active Galaxies in the Spitzer Infrared 
Nearby Galaxies Survey (SINGS). Thus the environment of these translucent clouds is more efficient 
at producing rotationally excited H2 per PAH-exciting photon than the disks of entire galaxies. Ex- 
citation analysis finds that the S(l) and S(2) emitting regions are warm (T > 300 K), but comprise 
no more than 2% of the gas mass. We find that UV photons cannot be the sole source of excitation 
in these regions and suggest mechanical heating via shocks or turbulent dissipation as the dominant 
cause of the emission. The clouds are located on the outskirts of the Scorpius-Centaurus OB associa- 
tion and may be dissipating recent bursts of mechanical energy input from supernova explosions. We 
suggest that pockets of warm gas in diffuse or translucent clouds, integrated over the disks of galaxies, 
may represent a major source of all non-active galaxy H2 emission. 

Subject headings: ISM: molecules— ISM: individual (DCld 300.2-16.9, MBM 33, LDN 1780)— ISM: 
lines and bands — infrared: ISM 



1. INTRODUCTION 

The kinetic temperature of all cold (T < 100 K) dif- 
fuse interstellar clouds is the result of a balance between 
heating and cooling processes. The H i and II2 receive 
kinetic energy from electrons ejected from dust grains 
and polycyclic aromatic hydrocarbons (PAHs) by inter- 
stellar far-ultraviolet (FUV; 6eV< hv < 13.6 eV) pho- 
tons, and transfer that energy to all other gas phase 
constituents. Clouds are cooled when inelastic colli- 
sions lead to spontaneous emission of spectral line radi- 
ation. Despite its being the most abundant constituent 
of molecular clouds, the role of H2 in directly cooling 
the diffuse IS M via its own line emission is uncertain. 
ISpitzed (fl94l was the first to estimate the significance 
of the pure-rotational {v = — 0) transitions of H2 
in cooling interstellar matter. Since that time, it has 
been realized that other species (C 11 in atomic gas and 
CO in UV-shielded molecular gas) were much more ef- 
ficient at converting kinetic energy to line emission, de- 
spite their much lower abundances (~ 10~*) relative to 
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hydrogen (" Wolfire. McKec. Hollenba ch. fc TielensI 120031 : 
iHollenbach. Takahashi. fc Tielensijgm 

Because of the small moment of inertia of H2, its rota- 
tional energy levels are difficult to excite in the cold ISM. 
The lowest-lying transitions are the electric quadrupole 
pure-rotational lines. Because H2 lacks an electric dipole 
moment, the J = 1 level is unable to radiate, so the first 
radiating level is J = 2. The first three rotational transi- 
tions are thus S(0) (J = 2-0) at 28.2/im, S(l) (J = 3-1) 
at 17.0/im, and S(2) (J = 4 - 2) at 12.3^m. All were ob- 
servable with the Spitzer InfraRed Spectrograph (IRS). 
The lowest transition S(0) has hu/k = 510 K. The S(2) 
transition occurs from a rotational level {J — A) that is 
about 1700 K above the ground state. 

Knowledge of the distribution of physical conditions 
along a given line of sight is required for the proper in- 
terpretation of any observation of H2 lines. This is es- 
pecially true for path lengths spanning many kpc such 
as the disks of galaxies. Since these quadrupole transi- 
tions require high temperatures to excite by collisions, 
the presence of H2 lines in the ISM of galaxies is usu- 
ally thought to be the result of either collisional exci- 
tation in warm H2 or UV pumping. Collisional exci- 
tation is expected to dominate in shock waves, while 
UV pumping dominates the excitation in photodisso- 
ciation regions (PDRsI (eg.. iShuU fc HoUenbacbl 119781: 
iBertoldi et all 120001 : iMever et al.ll2001[ ). In other cases, 
however, the source of H2 excitation is not so clear. In the 
recent Spitzer Infrared Nearby Galaxies Survey (SINGS; 
iKennicutt et al.l I2003D measu rements of H2 in normal 
galaxies. iRoussel et al.l (|2007[ ) concluded that FUV-rich 
photodissociation regions (PDRs) at the interfaces be- 
tween bright H 11 regions and molecular clouds were the 
most likely source of the emission. The authors did note 
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Table 1 

Observed H2 and Dust Surface Brightness 



H2 Pure-Rotational Surface Brightness^ 



Position 


RA (2000) 
(hh:mm:ss.s) 


DEC (2000) 
(dd:mm:ss) 


S(0) 


S(l) 


S(2) 
(10- 


S(3) 
■iOwm-2 


S(4) 
sr-i) 


S(5) 


S(6) 


(10-7) 


TIR'^ 
(10-6) 


DCld 300.2-16.9 (C) 


11:48:24.4 


-79:18:00 


36.74 


13.5 


12.8 


< 4.0'* 


< 6.8 


< 4.7 


< 8.0 


2.436 


1.64 






±0.16 


±0.4 


±0.5 










±0.007 


±0.16 


DCld 300.2-16.9 (B) 


11:52:08.3 


-79:09:33 


29.43 


13.5 


8.9 


< 3.8 


< 5.8 


< 5.3 


< 11.2 


2.745 


3.4 






±0.16 


±0.4 


±0.4 










±0.015 


±0.3 


DCld 300.2-16.9 (A) 


11:54:23.1 


-79:31:42 


24.03 


9.4 


7.7 


< 3.4 


< 4.0 


< 4.3 


< 5.4 


2.137 


1.51 






±0.15 


±0.4 


±0.5 










±0.018 


±0.15 


DCld 300.2-16.9 (D) 


11:55:33.8 


-79:20:54 


26.27 


9.7 


9.0 


< 4.2 


< 5.1 


< 5.0 


< 5.4 


1.714 


0.83 






±0.14 


±0.4 


±0.5 










±0.005 


± 0.08 


LDN 1780 (7) 


15:39:23.0 


-07:10:05 






6.9 


< 5.4 


< 3.2 


< 1.5 


< 3.2 


0.742 


1.47 










±0.6 










±0.003 


±0.14 


LDN 1780 (3) 


15:40:34.2 


-07:13:14 






20.9 


< 4.6 


< 5.4 


< 9.5 


< 7.6 


1.533 


2.7 










±0.6 










±0.002 


±0.3 



^ Integrated H2 surface brightness from spectrum fits, in units of 10^^*^ W m^^ sr^ ^ . 
^ ^7. 9, PAH = {^^u)7.9^im measured over the Spitzer IRAC 7.9/Am band. Units arc 10^' Wm^^sr 
Estimated total infrared surface brightness, in units of IQ- ^ Wm-^ sr^^. 



Upper limits are 3cr, as described in the text. 



that some of the H2 flux could originate in diffuse gas 
illuminated by the ambient FUV field (Gq ^ 1, in units 
of the average interstellar radiation field near the Sun) , 
which is in fact th e source of most of the IR emission in 
the same sample (jPraine et alJ [20071 ) . But their paper 
did not favor this interpretation of their data. 

A growing body of analysis finds that the H2 emission 
from the disks of galaxies does originate partly in diffuse 
material. Using FUSE measurements of FUV absorp- 
tion by H2 towards late B stars behind the Chamaeleon 
complex, iGrv et aD ()2002| ) found that the rotational lev- 
els could not be populated by only FUV photons, and 
required collisional excitation within anomalously warm 
regions. Indee d, further study of th e cirrus cloud in front 
of these stars (jNehme et al.ll2008af) has determined that 
the observed CH+ absorption also depends on the exis- 
tence of warm pockets of gas, since the formation of CH+ 
molecules needs an activation energy of AE/k — 4640 K. 
The first detection of pure-rotat ional H2 emission from 
diffuse Galactic regions was the iFalgarone et al.l (|2005[ ) 
ISO SWS measurement of the S(0) to S(3) transitions 
along a direction in the mid-plane of the Milky Way 
(MW: I, 6=26 ?46,0?09) As con cluded from the absorp- 
tion result of IGrv et all (|200a). the MW H2 emission 
turns out to be too bright to be caused only by FUV 
excitation. 

In this paper, we give additional evidence that low ex- 
tinction, low-UV environments can excite II2 molecules. 
In the following section we describe recent ultra-sensitive 
Spitzer Space Telescope observations of the 5-15/im (5~ 
38fim for four postions) spectrum towards translucent 
clouds at high Galactic latitude > 15deg) associated 
with infrared cirrus. In §3 we show detections of H2 S(0) 
and S(l) emission towards 4 positions and S(2) emission 
towards 6 positions. We estimate the physical conditions 
in the regions producing II2 line radiation. We compute 
ratios of H2 fiux to PAH feature strength and compare 
with the measurements of the same quantities made as 
part of the Spitzer Legacy SINGS survey of non-active 
galaxies. In §4 we review possible sources of excitation 
of the H2 molecules, and discuss the implications for the 



study of II2 emission from the disks of galaxies. In §5 we 
make our concluding remarks. 

2. OBSERVATIONS 
2.1. Spitzer Spectroscopic Measurements 

As part of a study of the heating and cooling of high 
Galactic latitude translucent clouds (HLCs), we spec- 
trally imaged 34 positions in six HLCs at 5.2-14. S^tm 
using the Spitzer Space Teles cope (iWerner et al.l [20041) 
Infrared Spectrograph (IRS) CHouc k et alJ 120041 ) . The 
data were taken under General Observer programs 20093 
and 40207. We detected six positions in two of the clouds, 
DCld 300.2-16.9 and LDN 1780 {I, b w 358?9, +36?9), in 
H2 S(2) pure-rotational emission at 12.3/im. The data 
for position A of cloud DCld 300.2-16.9 {Spitzer Request 
Key [SRK] 0014030848) were obtained on 2005 July 9; 
those for positions B, C, and D (SRK 0022715136) were 
obtained on 2007 July 24; and those for LDN 1780 posi- 
tions 3 (SRK 0022715648) and 7 (SRK 0022715392) were 
obtained on 2008 March 26. In addition, we imaged the 
14-38/im spectrum towards four positions in DCld 300.2- 
16.9 under Director's Discretionary Time program 491. 
These measurements (SRK 0028315392) were obtained 
on 2009 April 16, and yielded H2 S(0) 28.2^m and S(l) 
17.0yLtm detections towards all positions. The coordinates 
of the H2-detected positions are given in Table [TJ Those 
for the undetected positions are provided in Table [2] Al- 
though we present some of the dust measurements in this 
paper, a full analysis of the interstellar dust emission in 
the six-cloud sample is the subject of a future paper (In- 
gaUs et al, in preparation). We describe in this section 
the method by which we measured the S(0), S(l), and 
S(2) emission. 

Figure [1] shows the orientations and areas of coverage 
of the IRS subslits that we used to produce spectra for 
the four positions in DCld 300.2-16.9 and the two posi- 
tions in LDN 1780. The slit apertures are shown super- 
imposed on 100/Ltm emission maps from the IRAS Sky 
Survey Atlas (ISSA). The 5.2-14.5/im measurements of 
both clouds were made using the IRS short wavelength 
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Table 2 

Positions with H2 Upper Limits and Dust Surface Brightness 



H2 Pure-Rotational Upper Limits'' 



Position 


RA (2000) 


DEC (2000) 


S(2) 


S(3) 


S(4) 


S(5) 


S(6) 


r b 

J y g PA JJ 








(h.h:iiim:ss.s) 


(dd.:iiim:ss) 






Wm-2 






(i6-'7) 


MBM 12 (01) 




02:56:09.6 


+ 19:26:51 


39.3 


60.7 


90.4 


54.7 


85.6 


1.264±0.018 


MBM 12 (A 




509.1,-f636.4)'= 


02:56:19.9 


+ 19:26:04 


2.6 


4.1 


2.6 


8.4 


3.7 


0.573±0.007 


MBM 12 (A 




381.8,-1-509.1) 






2.2 


4.7 


2.5 


8.5 


5.8 


0.681±0.022 


MBM 12 (A 




254.6,-1-381.8) 






11.6 


40.2 


28.2 


15.6 


23.7 


1.09±0.03 


MBM 12 (A 




254.6,-1-381.8) 






11.6 


40.2 


28.2 


15.6 


23.7 


1.09±0.03 


MBM 12 (A 




127.3,-1-254.6) 






15.6 


78.0 


35.7 


23.6 


33.7 


0.88±0.05 


MBM 12 (A 


0,-1-127.3) 






12.4 


52.6 


24.9 


24.6 


23.6 


0.62±0.03 


MBM 12 (A 


-1-127.3,-254.6) 






14.0 


102.0 


32.3 


22.8 


28.8 


0.78±0.24 


MBM 12 (A 


-1-254.6,-381.8) 






9.0 


26.8 


17.5 


16.2 


18.6 


0.962±0.022 


MBM 12 (A 


-1-381.8,-509.1) 






2.9 


4.1 


3.8 


7.0 


3.6 


0.634±0.015 


MBM 12 (A 


-1-509.1,-636.4) 


" 


" 


2.5 


5.4 


2.7 


9.7 


4.9 


0.570±0.006 


MBM 12 (10) 




02:56:40.4 


+ 19:24:20 


24.2 


87.3 


72.4 


41.9 


57.6 


1.15±0.03 


MBM 12 (01 





ff) 


02:57:07.9 


+19:16:17 


14.8 


38.2 


38.0 


22.0 


29.9 


1.039±0.022 


MBM 12 (12) 


02:57:27.2 


+20:02:45 


33.7 


51.3 


79.9 


56.8 


61.7 


1.866±0.019 


MBM 12 (12 





ff) 


02:58:39.0 


+20:14:16 


45.0 


69.4 


106.7 


62.5 


85.0 


1.472±0.017 


MBM 28 (Af 




09:29:09.2 


+70:31:00 


35.5 


55.0 










DCld 300.2-16.9 (E) 


11:46:18.9 


-78:46:33 


1.3 


2.5 


3.1 


7.2 


5.6 


0.562±0.011 


LDN 183 (Wl 


-540,0) 


15:52:54.5 


-02:52:24 


2.7 


7.9 


6.7 


6.1 


7.6 


1.341±0.018 


LDN 183 (Wl 


-360,0) 






2.9 


6.6 


6.9 


3.8 


4.7 


1.527±0.015 


LDN 183 (Wl 


-180,0) 






3.8 


9.2 


10.0 


6.6 


7.6 


1.675±0.017 


LDN 183 (Wl 


0,0) 






3.2 


9.2 


11.2 


6.5 


5.0 


1.329±0.013 


LDN 183 (Wl 


-1-180,0) 






2.7 


4.0 


7.3 


13.9 


13.2 


1.033±0.009 


LDN 183 (Wl 


-(-360,0) 






2.7 


4.8 


6.9 


12.9 


13.2 


0.960±0.018 


LDN 183 (Wl 


-f 540,0) 






2.8 


5.0 


6.7 


13.3 


14.3 


0.902±0.009 


Stark 4 (SI; 


-480,0) 


23:45:07.3 


-71:42:47 


5.2 


9.1 


13.7 


12.3 


16.2 


< 0.340° 


Stark 4 (SI; 


-240,0) 






3.3 


5.4 


7.4 


6.1 


8.4 


0.113±0.016 


Stark 4 (SI; 0,0) 






2.1 


3.4 


5.3 


5.2 


6.4 


0.206±0.019 


Stark 4 (SI; -1-240,0) 






3.4 


4.7 


8.0 


5.3 


5.3 


0.178±0.006 


Stark 4 (SI; 


-1-480,0) 






7.7 


10.4 


15.5 


11.1 


14.0 


0.299±0.012 



^ Integrated H2 surface brightness upper limits (Scr) from spectrum fits, in units of 10 ^'^ Wm ^ sr ^. 
^ -f7.9,PAH = (i^Ii^)? .9fim mcasurod over the Spitzer IRAC 7.9/im band. Units arc 10^' Wm~^ sr""'". 
^ Offset from AOR center given as (ARA,ADEC), in arcseconds. 

^ The SL2 portion of the spectrum for MBM 28 (A) was corrupted by stray light and did not yield measurements for S(3) 
through S(6), nor for I?. 9, pah- 

^ Upper limit to /t.q.pah is 3 times the RMS in the band. 



MJy sr"^ 




Figure 1. IRAS Sky Survey Atlas images of the lOOfim emission tov^^ards (a) DCld 300.2—16.9 and (b) LDN 1780. Superposed are images 
of the SL (yellow) and LL {red) slit apertures used for our measurements. The direction tov^^ards the Lov^^er Centaurus Crux and Upper 
Scorpius OB associations are shown by arrows on the left and right panels, respectively. In these images, celestial North is up and Galactic 
coordinates are indicated by a dashed grid. 
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low resolution (SL) module. Each target was observed 
in standard staring mode for 8 cycles of 60 seconds each. 
Under staring mode, for each data collection cycle the 
SL2 subslit (5.2-7.7^m), and then the SLl subslit (7.4- 
14.5/im), are each placed on the target position at two 
offset "nod" positions, 1/3 and 2/3 along the long axis of 
the 57" X 3.7" subslit. When one subslit observed the tar- 
get, the other subslit observed a region of adjacent sky. 
The two subslits are separated by 22". [For a complete 
description of the SL geometry, see the IRS Instrument 
Handbook (SSC 2009).] We spent a total of about 32 
minutes imaging a 234" x 3.7" region around each target. 
Since each subslit observes a different region of adjacent 
sky, the area of overlap between the two subslits was only 
95" X 3.7". Nevertheless in what follows we combine all 
observations in the 234" window covered by either SLl 
or SL2 (yellow in Fig. [T]). The results do not change if 
we limited the measurements to within the 95" region of 
overlap, although the noise increases noticeably. 

The 14-38/im measurements towards 4 positions in 
DCld 300.2-16.9 were made using the IRS long wave- 
length low resolution (LL) module, whose slit covers a 
larger region of the sky than SL. Each of the four targets 
(plus the background position — see below) was observed 
in staring mode for 24 cycles of 30 seconds each. As for 
the SL observations, the short (LLl subslit; 14.0-21.3/xm) 
and long (LL2 subslit; 19.5-38.0/im) wavelength portions 
of the spectra were imaged on different portions of the 
sky. In the case of LL, however, we only combined the 
data for the two subslit positions that overlapped on the 
sky, yielding about 24 minutes per position and covering 
a 280" X 10.6" region (red on Fig. [T] a). To increase 
the likelihood of scheduling these measurements, we did 
not constrain the observation dates, and therefore the 
SL slit position angles were not replicated exactly by the 
LL measurements. This resulted in SL and LL effective 
apertures that were well-matched in size, but only over- 
lapped a relatively small portion of sky, near the centers 
of the imaged regions (see Fig. [Ij. In other words, we 
have chosen to minimize random noise (maximize inte- 
gration time per area of sky) at the expense of systematic 
error due to mismatched position angle, assuming that 
the emission is not a strong function of position. 

Our clouds all lie in regions of relatively bright Zo- 
diacal emission [according to the Spitzer Planning Ob- 
servations Tool (SP0Tj3 background estimator, the Zo- 
diacal foreground at 7.7/im ranges from 5-15 MJy sr~^]. 
To measure this sky foreground as well as correct for 
short term detector drifts, every observation (or con- 
current group of observations) of a given cloud was 
accompanied by an off-source "background" measure- 
ment with the same integration parameters. The back- 
ground position was taken to be a local minimum in a 
5° X 5° IRAS Sky Survey Atlas (ISSA) lOO^im map cen- 
tered on the cloud. The following sky background po- 
sitions were observed (equatorial J2000): DCld 300.2- 
16.9 (11:41:33.7,-78:11:37); LDN 1780 (15:33:50.9,- 
07:28:26); MBM 12 (02:51:36.0,-^19:10:03); MBM 28 
(09:38:34.9,+71:13:58); LDN 183 (15:48:48.8, -00:34:35); 



http:/ /ssc. spitzer. caltech.edu/warmmission/propkit/spot/ 



and Stark 4 (23:52:41.1,-71:12:45). 



2.2. Data Reduction 

Data processing started with the Basic Calibrated 
Data (bed. fits) spectral image files from the S18.7 
Spitzer IRS calibration pipeline. As mentioned above, 
to maximize the signal to noise ratio in the SL data, 
we combined all 32 spectral images towards a given 
position, covering a 234" x 3.7" region on the sky. 
We used an outlier-resistant computation of the mean 
(resistant jnean. pro from the IDL Astronomy User's 
LibrarjQ) of the 32 measurements for each pixel. For 
the LL data, we combined similarly the 48 spectral im- 
ages covering the inner 280" x 10.6" region covered by 
both subslits. To remove Zodiacal foreground emission 
and residual dark current, and to minimize the effects 
of "rogue" pixels that vary on timescales greater than 
a few hours, we subtracted a mean image of the nearby 
sky (see above) from each mean source image. We used 
IRSCLEAJ^ to interpolate over any remaining rogue 
pixels in each of the sky-subtracted mean images. We ex- 
tracted full slit extended source calibrated spectra from 
the resulting images using the Spitzer IRS Custom Ex- 
traction tool (SPICE) O 

We combined the surface brightness spectra for SL or- 
ders 1 and 2 into a single spectrum covering the wave- 
length range 5.79-14/Ltm. To do this we re-gridded the 
AA ~ 0.04/Ltm SL2 spectrum to match the 0.06^m SLl 
wavelength grid using linear interpolation, and averaged 
SLl and SL2 values where the two orders covered the 
same wavelengths. Before combining, the edge of each 
order was trimmed by from 2-15 spectral elements to 
eliminate regions of poor responsivity. This process was 
repeated for the LL data, giving a combined LL spec- 
trum with a resolution of 0.17/im covering the wavelength 
range 14.2-35.8//m. To obtain a 5.79-35.8/im spectrum 
for the four positions with LL measurements, the SL and 
LL spectra for these four positions were concatenated to- 
gether, retaining t he different resolut ions. We used the 
PAHFIT package (jSmith et al.l 120071 ) to fit simultane- 
ously a set of Polycyclic Aromatic Hydrocarbon (PAH) 
dust features, gas phase lines (including the H2 pure- 
rotational lines), dust continuum emission, and an ex- 
tinction curve. The fits were weighted by the error spec- 
trum produced by the Spitzer IRS pipeline and SPICE. 
Uncertainties for all fit parameters were computed by 
PAHFIT using full covariance matrices for the multicom- 
ponent fits. A proper fit to the data required modifica- 
tion of the defauh PAHFIT parameters. The 12.69^m 
PAH feature seems to be considerably wider in these 
spectra than expected from PAHFIT default parameters, 
resulting in spurious detections of the [Ne 11] 12.8/im line 
when the default narrow 12.69/im width was used. This 
left a "knee" of emission in the model spectrum that was 
not visible in the data. We found that increasing the 
fractional FWHM of the 12.69^m Drude profile, from 
7^(12.69) = 0.013 to 0.042 (the same as for the 12.62/im 

* http://idlastro.gsfc.nasa.gov/ 

^ http://ssc.spitzer.caltech.edu/dataanalysistools/tools/irsclean/ 
http: / / ssc.spitzer.caltech.edu / dataanalysistools / tools / spice / 
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Figure 2. Spitzer IRS observed spectra (filled circles with error bars) for 4 positions in DCld 300.2-16.9 in front of tiie Cliamaeleon complex. Superimposed on the spectra are fits 
(solid black line) made using the PAHFIT program (Smith et al. 2007), that are each the sum of three dust continuum components (brown), a 5000 K stellar blackbody (green), gas 
phase spectral lines (red), and a set of PAH features (blue), and account for silicate extinction. The thick gray curves indicate the total continuum emission. The inset figures show 
the residuals in the data for a portion of the spectrum (solid histogram), after subtracting off the fits to the dust continuum and PAH features and correcting for extinction. The 
summed fits to all gas phase spectral lines are superimposed in red on the inset plots. All inset plot windows cover the same range in wavelength as the region of the main plot window 
they occupy. The size of each inset window in its main plot window is given by a gray rectangle overlaid on the full spectrum. 
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feature), removed the spurious knee and eliminated the 
[Ne ii] "detection." 

For the S(3) through S(6) H2 hues, which were not de- 
tected, we estimated 3ct upper hmits to the integrated in- 
tensity as follows. First, we subtracted the fit to all PAH 
features and dust continuum components from the ob- 
served spectrum. We divided the result by the PAHFIT- 
derived extinction curve (which in most cases did not 
affect the S(2) or 7.9/zm intensities by more than 10%). 
For a given H2 line we then determined armsi the root 
mean square of this residual extinction-corrected spec- 
trum in a wavelength region centered on the line center 
wavelength, Aq ± 4 x FWHMa . Since all molecular lines 
are unresolved by the IRS, the expected full width at half 
maximum (FWHMa) of the line is determined only by 
instrumental parameters fsee lSmith et al.l[2007L for more 
details on the H2 line parameters). For a gaussian line 
profile with peak intensity Ix,peak the integrated intensity 
is given by 



/ = / hd\ - 



27r(FWHMA) 



X,peak 



2V2hr2 

1.064 (FWHMA)/A,peafe. (1) 



Thus the 3cr integrated surface brightness limit is deter- 
mined by substituting 3 arms for Ix.peak in Equation [T] 



I fi / Ao \ /FWHMa 
w -2 -1 < 3-191 X 10-g ^ 

3 (Jrms 



MJy sr^i 



(2) 



2.3. Dust Emission 



To constrain the amount of interstellar radiation pro- 
cessed by cirrus clouds, we synthesized broadband mea- 
surements of the dust emission. For all positions, we 
computed /r.g.PAH = (j^^iy)7.9.PAH, the PAH surface bright- 
ness in Spitzer Infrared Array Camera (IRAC) band 4, 
which is centered on 7.9/^m. To do this we integrated the 
model PAH spectrum (given by the sum of fits to PAH 
features in Figures [5] and El weighted by the IRAC 4 filter 
resp onse, and m ultiplied the result by the color correc- 
tion (|SSCll2010al ) to give the specific intensity ( MJy sr'^) 
that would be measured photometrically on an IRAC im- 
age. Multiplying by the effective frequency of the 7.9/xm 
array gives the surface brightness (Wm~^sr~^). This 
procedure was implemented using the Spitzer Synthetic 
Photometry software n Error bars for /t.o^pah are de- 
rived from the covariance matrix for the fit PAH pa- 
rameters. For the position in Table [2] that did not have 
detectable PAH features we computed upper limits as 

3 times the RMS spectrum integrated over IRAC band 

4 (the spectrum was squared, integrated over the band, 
and multiplied by the color correction; the square root 
of this was multiplied by the effective frequency) . 

We made a crude estimate of the total infrared (TIR) 
surface brightness to compare with the total H2 emission 
from the S(0) to S(2) transitio ns. We conve r ted th e for- 
mula given in Equation (22) of iDraine fc Lil (|2007f ) from 

http: / / ssc.spitzer.caltech.edu/dataanalysistools / cookbook /lO/ 



spectral luminosity to surface brightness: 

TIR « 0.95 {iyl^)7.9 + 1.15 {1^^)24 

+ {vIv)7i + {vIu)wo- (3) 

The V values are the effective frequencies (in Hz) of the 
Spitzer IRAC 7.9 and MIPS 24, 71, and leO^im band- 
passes, and In is the specific intensity measured with 
these instruments. 

The surface brightness {vlv)7.9 was obtained as above. 
For the four positions with LL spectra in DCld 300.2- 
16.9, {vIi,)2A was derived similarly using the LL spectrum 
and the MIPS 24/im filter respon se, times an appropri- 
ate color correction (|SSO l2010bD . For the LDN 1780 
positions, which did not have LL spectra, we estimated 
the 24^m contribution to TIR as follows. Taking the 
measured average value of (t/Ji^)7.fl/ TIR = 0.24 ± 0.03 
for our sample, we used Figure 15 of IDraine fc O (l2007h 
(second panel) to infer a dust mass fraction in PAHs of 
9pah = 4.6%. For this value of ^pahi the third panel of 
the Figure predicts (i//y)24/TIR w 0.06. This assumes 
band emission relative to TIR does not vary as a function 
of illuminating radiation field strength U (i.e., U < 10, 
which probably holds for all cirrus clouds). 

Lacking 71 and 160^m observations towards the DCld 
300.2-16.9 positions, we used IRAS Sky Survey Atlas 
(ISSA) 60 and lOO/um measurements, respectively, as 
proxies for the missing MIPS photometry. One of the 
HLCs in the fuh six-cloud sample, MBM-12 (LDN 1457) 
(undetected in H2 S(2) emission and not observ e d wit h 
LL), was mapped with MIPS bv lMagnani. et all (|2006[) . 
We computed the following rough scaling relations: 



/,,7i = (2.30 ± 0.08) /,,6o 
/.,i60 = (4.98 ± 0.13) /,,ioo. 



(4) 
(5) 



These were derived from robust linear fits to scatter plots 
of vs. /,y,60 and Iumo vs. /i.,ioo for MBM-12, using 
images in the given wavebands resamplcd to the same 
pixel grids. The errors in TIR are dominated by the 10% 
accuracy of Equation[3l as determined from model clouds 
heated by starlight with intensities between 0.1 and 10 
times the average interstellar value (|Draine "¥Lil[200l . 

3. RESULTS 

3.1. Measured Surfaee Brightnesses 

The complete 5.8-35.8/im SL + LL spectra for the 
DCld 300.2-16.9 positions are displayed in FigureEl The 
5.8-14.0^m SL spectra for the LDN 1780 positions are 
shown in Figure |31 We overlay on the spectra the fits to 
the grain and gas phase emission, and label the wave- 
lengths of all lines in the model (regardless of detec- 
tion). We detected emission in the (■i;=0-0) 28.2/im S(0) 
and 17.0/im S(l) transitions of H2 towards all four DCld 
300.2-16.9 positions observed with the Spitzer IRS Long 
Low (LL) spectrograph. We detected emission in the 
12.3/im S(2) line at the > 3a level in only six of the 34 
HLC positions observed with the Short Low (SL) spec- 
trograph: the four in DCld 300.2-16.9 observed with LL 
and two in LDN 1780. The 3a detection threshold for 
S(2) was between 0.7 and 5.4 x 10"^° Wm"^ sr"^ for the 
34 observations. For the cloud positions with detectable 
S(2) emission, the S(0), S(l), and S(2) surface bright- 
nesses and their uncertainties (computed by PAHFIT) 
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Figure 3. Same as Figure[2] except that these spectra are for 2 positions in LDN 1780, and only cover the SL wavelength range 5.8-14.0/^m. 
Inset boxes cover a diiTerent range in wavelength and surface brightness per wavelength than in Figure [2] (but cover the same ranges in 
these two panels). 



are listed in Tabled! columns 4 through 6. 

The S(3), S(4), S(5), and S(6) lines of H2 were all too 
faint to be detected: the integrated line strengths de- 
rived from attempting to fit the lines at the expected 
wavelengths were smaller than the corresponding errors 
in the fits. Upper limits (3cr) to the surface brightness 
in these transitions (Equation [2]) average about 4.2, 5.0, 
5.1, and 6.8 x 10~^° Wm~^ sr~^, respectively. We print 
the actual limits in columns 7-10 of Table [TJ For the 
positions without detectable H2 emission, upper limits 
to the S(2) through S(6) intensities are listed in columns 
4-8 of Tabled 

We list in Table [T] the dust emission surface bright- 
ness: /y.g^pAH Is lu columu 11, aud the TIR estimates are 
in column 12. For the positions without detectable H2 
emission, column 9 of Table [2] provides /7.9.PAH measure- 
ments. 

3.2. Molecular Hydrogen Excitation 

We can use an excitation diagram to estimate the phys- 
ical conditions in the regions producing H2 rotational 
line emission in high-latitude translucent cirrus clouds. 
In Figure |4l we plot for all the Table [T] measurements 
the column density in rotational upper level J divided 



by the statistical weight, Nj/gj, as a function of Ej/k, 
the upper level energy expressed as a temperature. The 
column densities are derived from the surface brightness 
assuming optically thin line radiation: 



N.J 



A 



4:7^1 J, J -2 



J,J-2 



7-2 



7, ,7-2 



(6) 



is the 



Here, Ijj-2 is the surface brightness and A 
Einstein spontaneous emission coefficient in the transi- 
tion from upper state J to lower state J — 2. The sta- 
tistical weights for rotational level J are gj_^^^„ = 2J +1 
for even J (para-H2) and 5,7, odd = 3(2J + 1) for odd J 
(ortho-H2). In what follows we assume that excitation 

of vibration al levels is negligibl e. 

Following iDraine fc Bertoldil (|1996[ ) , we approximate 
the II2 rotational populations by a thermal distribution. 



Nj/gj = a eM-Ej/kTr), 



(7) 



where Tr is the "rotational temperature" describing the 
populations. Under conditions of local thermodynamic 
equilibrium (LTE) Equation[7]becomes a Boltzmann dis- 
tribution, and a — N (112) / Z (Tr) , the ratio of total H2 
column density to the partition function. Under non- 
equilibrium conditions, the ortho- and para-Il2 can take 
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Figure 4. Excitation diagram of H2 pure-rotational emission towards high-latitude cirrus clouds. Circles represent the observed measure- 
ments of DCld 300.2-16.9 (red) and LDN 1780 (blue). Transitions S(3) and higher all have Scr upper limits, shown as downward arrows. 
The abscissae for each point have been shifted slightly for clarity. Superposed on the diagram are seven cloud models: (1) a thermal 
distribution (Twarm = 300 K) with parametrized ortho/para ratio 7 = 0.70 (gray); (2) a thermal distribution (Tcoid = 85 K) with 7 = 1.72 
(gray); (3) emission from a PDR illuminated on both sides by a UV flux equal to the average local interstellar value (G?o = 1: IHabin3 
I196SI 1. volume density Uh = 30 cm""^, and visual extinction Av = 0.3 mag (dotted); (4) PDR emission with Go = 10, Uh = 100 cm~^, 
and Av = 2.0mag (dashed); (5) emission from 8 MHD shocks traveling at Vs = 8 kms~^ in gas with Uh = 50 cm""^ and magnetic field 
B = 7/iG (dot-dashed); (6) emission from 1333 magnetized vortices with ion-neutral drift velocity = 3.5 kms"^, Uh = 30cm~'^, and 
B = 10 ^iG (dot- dot- dashed); (7) the sum of the Go = 1 PDR and the vortex model (orange); (8) the sum of the 85 K and 300 K thermal 
distributions (solid black); and (9) the sum of a PDR w i th Gp = 0.6 plus a warm component with T = 391 K, with non-LTE excitation. 
Models (3)-(6) are taken from Fig. 3 of IFalgarone et ahl I I2005I1 . Model (9) is taken from Draine & Ingalls (in preparation). 



on separate a values, aortho and a^^,^. One way to char- 
acteri ze the relative population s in ortho and para states 
is the iDraine fc Bertoldil (|1996f) ortho to para ratio 



7 = 3a„tho/apa 



(8) 



As those authors note, 7 is not necessarily the col- 
umn density ratio of ortho- and para-H2, but rather a 
parametrization of that ratio that includes the ratio of 
partition functions of the two types of molecular hydro- 
gen. It has the advantage that under conditions of LTE 
(floitho = Opara), 7 = 3, rcgardlcss of temperature. 

The observed excitation diagram can be used to es- 
timate physical conditions in the S(0) to S(2)-emitting 
gas. Using Equation [71 we note that for constant val- 
ues of Tr and a, a plot of Nj/gj versus Ej/k with 
logarithmically scaled ordinate will yield a straight line 
with slope equal to — l/T^ (multiplied by l/logio(e) if 
the ordinate is log base 10). As is typical for other 
measurements of H2 rotational populations in diffuse 
gas (eg., f Rachford rgTall 120011: iFa garone efall 120051 : 
iGoldsmith et al.. .20101 : iJensen et al.ll2010[ ). our data do 
not fall on a single straight line. Here we simply de- 
termine two excitation temperatures consistent with the 
three lines available. The result is somewhat arbitrary 



and not unique, but nevertheless is instructive. 

We find that a "cold" component with excitation tem- 
perature T^oid = 85 K and a "warm" component with 
T'warm = 300 K arc consistent with the data. This simple 
model, with N(H2) = 2.6 x 10^^ cm~^ in the cold com- 
ponent and iV(H2) = 3.1 x 10^*cm~^ in the warm com- 
ponent, reproduces the level populations in the J — 2,3, 
and 4 levels (corresponding to the S(0), S(l), and S(2) 
transitions) and the upper limits in J = 5 and above. 
We overlay on Figure [l] the excitation diagram for each 
thermal component (gray curves), as well as the sum of 
the two components (solid black). The model predicts 
a J = 5 population that is right at the upper limit we 
have measured. This is not a concern, because the criti- 
cal densities for J — 3, A, and 5 are high enough that the 
J = 5 level will be sub-thermally excited (at T = 70K, 
n„,.(H2) = 900, 2 X 10"^, and 3 x 10^ cm'^ for J = 3, 4, 
and 5). Thus the gas temperature would actually have to 
be higher than T = 300 K to reproduce the populations 
of J = 5 predicted by our thermal model. 

In the model, we assume an ortho to para column den- 
sity ratio of iVortho/A^para — 0.7 for both the cold and 
warm components. This is based on the observed col- 
umn density ratio between the J — I and levels to- 
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Figure 5. Ratio of tlie H2 S(2) to 7.9/xm PAH flux as a function of the Ld/im surface brightness, for Milky Way and extragalactic sources. 
Red circles (A-D) a re our DCld 300.2-16.9 HLC cirrus detections; blue circles (3,7) are the LDN 1780 detections. The diamonds are 
IRoussel et al.l 1 120071 ) data for SINGS galaxies with non-active nuclei (filled diamonds are dwarf galaxies and open diamonds are galaxies 
with purely star-forming nuclei). The green circle (MW) is an estimate for thelRilgaronc et al. (2005) Milky Way galactic plane line of sight 
{Av = 30). Dashed ellipses (pinched due to the logarithmic scaling of the axes) indicate the average ±l<j for two samples of high-latitude 
clouds: (1) the 6 S(2)-detected positions; and (2) the full set of 34 positions observed with SL, including non-detections. 



wards HD 102065 ()Nehme et al.ll2008H) . about 3?4 away 
from the nearest position in Table [1] It seems reason- 
able to take the same ortho/para column density ratio 
for both components, on the theory that it is the cool 
gas that occasionally gets heated to T = 300 K, keep- 
ing the relative amounts of ortho and para H2 constant. 
The parametrized ratio, 7, which measures departures 
from equilibrium (defined by 7 = 3), is obtained by sum- 
ming the predicted odd and even J Boltzmann factors 
over all possible rotational levels (including J — and 1, 
which do not radiate), and taking the ratio. This yields 
7coid = 1-72 and 7warm = 0.70 for our thermal model. 

The two-temperature model (T^om = 85 K and T„.,„„ = 
300 K) presented here is an over-simplification. In reality 
there will be a range of temperatures present, with much 
of the H2 column density at T < 85 K (based on FUSE 
absorption measurements towards 38 translucent lines of 
sight, the mean J = 1 to J = rotational temperature 
is 67 K; see iRachford. et"all [20091 ) . and some material at 
temperatures between 85 K and T„^^^ required to account 
for the observed excitation of J — 2. 

3.3. Comparison with Measurements of the Disks of 
Galaxies 

To answer the question of whether translucent gas is 
a significant source of galaxian II2 emission, we com- 
pare the H2 and du st emission measurements with the 
IRoussel et al.l (|2007D SINGS galaxy sample and a mea- 
surement towards the mid-plane of the Milky Way. We 



plot the flux ratios F(S2) /-F^.q^pah as a function of the sur- 
face brightness /t sj^pah for the 6 HLC positions in Figure 
m (Since the solid angles of the H2 and PAH measure- 
ments are the same, the surface brightness ratio equals 
the flux ratio.) 

We also plot the SINGS extragalactic measurements 
of F(S 2)/f,. 

9.PAH -^7. 9. P AH iri galaxies without active 
nuclei (|Roussel et al.ll2007l ). We derived 7.9fj,m surface 
brightnesses for this sample by dividing the reported 
Rouss eFet al.l IRAC band 4 fluxes corrected for stellar 
emission by the solid angle for each galaxy observation. 

We add a final data point to Figure [HI the 
ISO SWS detection of S(2) along a direction in the 
mid-plane of the Mil ky Way (MW; 6=26?46,0?09; 
iFalgarone et all |2005| ). Lacking an IRAC 7.9^m 
measurement or its equivalent for this position, we 
scaled the IRAS 12p,m measuremen ts (IRIS processing; 
iMiville-Deschenes fc Lagachd l2005l ) for the MW posi- 
tion to the average 7.9/j,m/12/j,m ratio derived from HLC 
synthetic photometry (7.9/^m) and IRIS m easurements 
(12um). We fit a least squares bisector line (jlsobe et al.l 
Il990f ) to the 7.9 and 12/im measurements and used the 
slope as the nominal scale factor, /7.9,PAH/-fi2,iRis7 to mul- 
tiply by the 12/im MW measurement. We estimated the 
error in the predicted 9 pah(MW) value using the stan- 
dard deviation in the six measured HLC values of the 
-^7.9,pAH/-fi2,mis ratio (after subtracting the bisector fit y- 
intercepts). Multiplying this by /i2,iris(MW) gives an 
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Figure 6. Ratio of the total rotational H2 flux to TIR flux, as a function of the TIR surface brightness. Red circles (A-D) are our DCld 
300.2-16.9 cirrus detections, blue circles (3 and 7) are lower limits for LDN 1780, and diamonds arc Rousscl et al. (2007) data for SINGS 
galaxies with non-active nuclei (see Fig. [Sjl. A power law fit to the DCld 300.2-16.9 points (gray line) yields F(SO + SI + S2)/Ftir = 
2.0 X 10-*TIR-0-9. 



error in /t.o pah (MW). This error is probably an under- 
estimate, since it ignores possible systematic differences 
in the PAH spectra of HLC and MW populations. 

The six HLC positions have higher values of the 
F(S2)/F7 9 PAH ratio than the Milky Way mid-plane and 
all but three of the SINGS non-active galaxies. The 
average F(S2)/i^7 g pAu rati o for the dete c ted sa mple is 
(6.7 ± 1.6) X 10-3.' For the lRoussel etall (l2007l) SINGS 
sample the ratio averages (2.9 ± 0.6) x 10"^ for dwarf 
galaxies, (1.44 ± 0.09) x 10"^ for galaxies with purely 
star-forming nuclei, and (1.8 zt 0.2) x 10"^ for the com- 
bined set. For the iFalgarone et al.l ()2005[ ) Milky Way 
position, the ratio is (1.0 ±0.6) x 10"^. For the complete 
34-position HLC sample, 27 lines of sight showed PAH 
emission in the 7.9/xm band but no S(2) emission (Table 
[2]). If we include all 7.9/um detections in the HLC aver- 
age, the ensemble has a ratio of (2.0±0.8) x 10"'^, which is 
identical within the error bars to the full SINGS average 
and the Milky Way data point. (We estimated the aver- 
age S(2) intensity for the complete sample using a spec- 
trum produced by coadding the background-subtracted 
IRS spectral images for all 34 positions.) In the following 
section we discuss the possibility that translucent cirrus 
comprises a significant component of galaxian H2 rota- 
tional emission. 

A useful diagnostic of the relative strength of H2 emis- 
sion to absorbed interstellar radiation is the ratio of the 
flux in the S(0), S(l), and S(2) lines to the total infrared 
flux, F(SO + SI + S2)/i^Tm- We plot in Figure [5] this ra- 



tio as a function of /tir for the four DCld 300.2-16.9 
positions with LL spectra, as well as the two lower limits 
from the LDN 1780 S(2) data. Here the DCld 300.2-16.9 
points lie well above the extragalactic measurements. 
The average value of F(SO -t- SI -I- S2)/i^Tm for DCld 
300.2-16.9 is (3. 5 ±0.9) x lO'^. For the SINGS data, the 
ratio is (3.7±0.6) x lO"" for dwarf nuclei, (4.8±0.2) x 10"'' 
for pure star-forming nuclei, and (4.5 ± 0.8) x lO"** for 
the complete sample. It is notable that for the LDN 1780 
positions, F(S2)/Ftir alone is of the same order of mag- 
nitude as i^(SO + SI + S2)/FTm for the SINGS galaxies. 
Unlike extragalactic H2 emission, the H2 emission for 
HLCs seems only loosely correlated (if at all) with the 
infrared flux: a power law fit to the DCld 300.2-16.9 data 
in Figure E] has 1(30 + SI + S2) cx . 

4. DISCUSSION 

We have detected molecular hydrogen rotational emis- 
sion towards six of 34 translucent high Galactic latitude 
cirrus positions. The (u = - 0) S(0), S(l) , and S(2) 
transitions of H2 originate from energy levels that are 
510 K, 1015 K, and 1682 K above the ground state. How 
are these levels excited in the translucent ISM, which is 
embedd ed in cold H i gas with a median temperature of 
- 70 K (|Heiles fc Troland 2003)? Why do some locations 
show S(2) emission and others do not? Do we expect 
a significant amount of H2 cirrus emission throughout 
the Galaxy, and by extension in other normal galaxies? 
These are the main questions we will try to address here. 
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4.1. The Environment of the Clouds 

We first examine the environment of the two clouds we 
have detected, DCld 300.2-16.9 and LDN 1780, focusing 
on aspects of the clouds that may explain the H2 emis- 
sion we have observed. Both clouds are thought to be 
part of the cold neutral medium (CNM) that defines the 
boundary of the 30-100 pc diameter Local Cavity (LC) 
of lo w-density (n < 0.1 cm~'^) gas in which the S un re- 
sides (|^iih]lLan^e^Z^rgSyISRaimon^ The 
clouds also both happen to be on the interface between 
the LC and the superbubble surroun ding the Scorp ius- 
Centaurus (Sco-Cen) OB association (jde Geuslll992[ ). al- 
beit on opposite sides. Southeast of the superbubble 
center (in Galactic coordinates), DCld 300.2-16.9 is lo- 
cated along the line of sight to the Chamaeleon molecular 
cloud, but is probably closer than the ^150 pc distance 
of the main complex (iMizuno et al.ll200ll iNehme et al.l 
I2008aj) . IMizuno et al.l (12001') deduced that the cloud is 
at a distance of 70 ± 15 pc from the Sun, based on the 
location of T Tauri star T Cha, towards which the cloud 
peak of CO emission has the same radial velocity. The 
cloud is thus ^ 30 pc from the center of the Lower Cen- 
taurus Crux (LCC) subgroup of Sco-Cen. Northwest of 
the superbubble center, LD N 1780 (£, b = 358?9, 36?9) is 
about 100 pc from the Sun ()Francdll989HLallement et al.l 
[2003) and is thus « 70 pc away from the ce nter of the Up- 
per Scorpius (US) s ubgroup of Sco-Cen (|Laureiis et al.l 
119951: Ide GeusI I1992D . Arrows pointing from the clouds 
to the centers of the LCC and US subgroups are drawn 
on Figure [1] Projected on the plane of the sky, LDN 
1780 is « 75° away from DCld 300.2-16.9, so the clouds 
are separated by at least 100 pc. 

In terms of their processing of the interstellar radiation 
field (ISRF), both clouds have 60/100/Ltm colors and C 11 
cooling intensities that are typical of translucent Galactic 
cirrus, and have been modeled as objects with visual ex- 
tinctions ~ 1 magnitude an d ISRF intensi ties Gg ^ 1 , 
in units of the iHabind (119681) fl ux (Jngalls et all 120021 : 
lJuvela. Padoan. fc Jime nez 2003| INehme et al.ll2008bD . 

We estimate that the Table [T] positions have extinc- 
tions of from 1 - 3 mag. For DCld 300.2-16.9, the 
IMizuno et all dMH) ^^CO map has integrated inten- 
sities Ico ~ lOKkms"^ near our positions, which 
yields a molecular column density of A^(H2) ^ 2 x 
10^^ cm~^, if we use a typic al conversion factor betwee n 
Ico and iV(H2) (eg., see iLiszt. Petv. k Lucas! I2010D . 
Using the well-established relationships Nh/Eb-v = 
5.8 X lO^^cm"^ and A^/Eb-v = 3.1 in diffuse ga s 
(IBohlin. Savage, fc Drakd [l97l: IRachford. et"all l2Q09l) . 
this gives « 1.1 mag (this is probably an underes- 
timate since it does not count H i). The LDN 1780 posi- 
tions have measured visual extinctions of A^ « 2 — 3 mag 
([del Burgo & Cambrcsy 2006,), which corresponds to a 
total gas column density of Nh ~ (3.7 — 5.6) x lO^-'^cm"^. 
Alternatively, the average 100/xm surface brightness to- 
wards all detected DCld 300.2-1 6.9 and LDN 1780 posi - 
tions is 13.4MJysr~^. Using the iBoulanger et afl (|1996D 
conversion A^(H i) ss 2 x 10^°cm~^ (7ioo/MJysr~^) 
yields a column density of 2.7 x 10^^ cm^^, or A^ = 
1.4 mag. 

The portion of DCld 300.2-16.9 in front of HD 102065, 
about 3?4 (4.2 pc) away from the nearest position in Ta- 



ble m has been studi ed in abso r ption of H2, CH+, and 
many o ther species bv lGrv et al.l ()2002D and lNehme et al.l 
([20083) • ^ model of the direction towards this star can 
reproduce most of the observa tions if A^ ~ 0.7 and 
Go 0.7 (INehme et al.ll2008 bD. However, the popula- 
tions in the H2 J > 2 rotational levels require excita- 
tion beyond that produced by any reasonable incident 
radiation field to explain the observed absorption. The 
authors concluded that a warm (T ^ 250 K), out-of- 
equilibrium ga,s is al so present along the line of sight. 
INehme et al.l (|2008a[ ) posited that the cloud's current 
state is the result of an interaction between a supernova 
blast wave from one or more stars in the LCC subgroup 
and the ambient ISM, about 2 — 3 x 10^ yr ago. 

LDN 1780 is also thought to have been affected pro- 
foundly by its proximity to the Sco-Cen association. The 
cloud is located between two near-parallel arcs of H i 
emission. The outer arc is assumed to be H i swept 
up by supernova explosions in the Upper Centaurus Lu- 
pus (UCL) subgroup of Sco-Cen, whereas the inner arc 
is believed to have been produced by the supernova e 
in the Upper Sco rpius (US) subgroup (|de Geu^[T99l ). 
iToth et all (|1995[ ) proposed a scenario in which a shock 
front propagating outward at a speed of 10 — 15 kms~^ 
from UCL has passed through LDN 1780. Since the loop 
is currently about 7pc away from the cloud, we estimate 
that the cloud's interaction with the shock probably oc- 
curred '-^ 6 X 10^ yr ago. Toth et al. noted that the H i 
emission from LDN 1780 peaks on the lower latitude dif- 
fuse side of the cloud facing the US subgroup and the 
stars of the Galactic mid-plane; the CO emission peaks 
on the dense side of the cloud further from these sources 
of radiation. They argued that the cometary shape of 
the cloud is due to the passage of the UCL shock wave, 
whereas the segregation between molecular and atomic 
material is due to an asymmetric UV radiation field origi- 
nating from the US s ubgroup and the Galactic mid-plane. 
iLaureiis et al.l (|1995[ ) estimated that Go might be as high 
as 3 on the diffuse side of the cloud, but may be only half 
as strong on the dense side. The H2 S(2) emission we de- 
tect is stronger towards position 3 (Table [Ij , which is 
close to the H i peak. We argue in §4.2 that the UV 
field is insufficient to produce the emission, and suggest 
instead that the S(2) emission results from the cloud's 
interaction with the supernova blast wave. 

LDN 1780 is unusual among molecular clouds in that 
it has a significant Ha surf ace brightness, correlated with 
tracers of column d ensitv (Idel Burgo &: Cambresvll2006l : 
iWitt et all IMS). iMattila. Juvela. fc LehtinenI (|2007f l 
suggested that this could be explained by scattering of 
backgrou nd Ha pho t ons b y dust grains, and the recent 
study bv lWitt et al.l (|2010[ ) strongly supports this inter- 
pretation. The C Oph and 6 Sco H 11 regions appear to 
be within a few tens of pc of the southern surface of LDN 
1780, and represent likely sources of background Ha. 

4.2. Sources of Molecular Hydrogen Excitation 

In interstellar cirrus, H2 may be excited into rotational 
emission via UV pumping by starlight; or via collisional 
excitation resulting from mechanical ener gy dissipation 
in magnetohydrodynam ic (MHD) shocks ((Draine et al.l 
Il983t iFlower et al.lll986l). coheren t vortices or intense ve- 
locity shears ('joulain. etaLllIggl iFalgarone et al.l [20051: 
iGodard et al...2009, ). or by interaction with cosmic rays 
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(jFerland et al.|[2008l) . All of these sources may be present 
at some level in DCld 300.2-16.9 and LDN 1780. 

4.2.1. UV Pumping 

The interstellar radiation field (ISRF) is the most ob- 
vious source of H2 excitation in cirrus clouds. FUV ra- 
diation heats the gas and dust and drive s the c hemistry 
in all translucent regions. Rou ssel et al.l (|2007[ ) found a 
correlation between H2 emission and infrared emission 
in the disks of SINGS normal galaxies. Since IR radia- 
tion traces the processing of starlight by dust, SINGS H2 
emission appears to be tied to the FUV radiation field. 
But for cirrus clouds, it is difficult to reproduce the ob- 
served II2 emission with only the ISRF. We overlay on 
Fig. m the predicted excit ation in two photodi ssocia- 
tion region (PDR) models (iFalg arone et al.H2005| ). The 
first model is immersed in the Habing (1968|) interstel- 
lar radiation field (Gq = 1) and has gas volume den- 
sity n = 30cm~'^, and visual extinction — 0.3 mag. 
The second PDR model has a larger UV flux, and is 
denser and thicker: Go = 10, n = 100 cm~^, and 
Av — 10 mag. The translucent HLCs in our study are 
probably somewhere in between these models. Both 
models fall short of reproducing our observed measure- 
ments, by at least a fa ctor of 10. In another H? fluo - 
rescence model, that of iBlack and van DishoedJ ()1987[ ). 
the S(2) emission we observe can only be reproduced 
if the radiation field was increased to Go ~ 30. This 
would, however, result in S(3) emission at a level we 
would have detected. Furthermore, there is no evidence 
that these HL Cs are heated by radiati on fields greater 
than Go » 3 (iLaureiis et a l. 1995; Ing alls et all 120021: 
lJuvela. Padoan. fc Jimened 12003; .Nehme et al.ll2008bD . 
So another mechanism must dominate the H2 pure- 
rotational excitation. 

We note that by far the brightest S(2) emission we 
have detected is towards position 3 of LDN 1780. This 
position has 3 times as much S(2) radiation as position 
7, and also has twice as much 7.9/im PAH emission as 
position 7, consistent with the PAHs on the southern side 
of the cloud (in Galactic coordinates) bei ng exposed to 
higher starlight intensities, as predicted bv lLaureiis et al.l 
([1995) . The ISRF would have to be increased by a factor 
of more than 100, however, to account for the H2 S(2) 
emission towards this position, which would boost the 
PAH emission by a similar factor (provided the density 
remained the same) and would result in detectable S(3) 
emission and much brighter S(0) and S(l) lines than we 
observe. The highest IRAS 60/100/im surfa ce brightness 
ratio measured towards LDN 17 80 is 0.275 (llngalls et al.l 
[200^ . which implies Go ~ 2-4. iWitt et aLlpOlOD found 
that for the entire cloud, grain temperatures are in the 
range 14.5-16.8 K, consistent with Go ^ 1. Therefore it 
is likely that most of the S(2) radiation represents the 
dissipation of mechanical energy that was concentrated 
on the south side of the cloud (like the passage of a su- 
pernova blast wave). 

4.2.2. Collisional Excitation 

H2 rotational transitions can be coll isionally ex- 
cited into emission via shock h eating (IDraine et al.l 
19831) or turbulent dissipation (jGodard et al.l 120091: 



What kind of shocks would be capable of powering the 
observed emission? The observed H2 S(0) to S(2) line 
intensities sum to 7(0-2) = I[S{0)]+ I[S{1)]+ I[S{2)] « 
5 X lO-^Wm-^sr-i (on average) at the 4 DCld 300.2- 
16.9 positions where emission is detected. If a fraction 
/(O — 2) of the shock power is radiated in these lines then, 
assuming the shock front fills the beam. 



7(0 - 2) 
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87r cos 9 



/(O - 2)pov'^ 
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for a strong shock, where 6 is the angle between the shock 
normal and the line of sight, po is the density of the 
pre-shock gas, and Vg is the shock propagation velocity. 
Thus, 
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Therefore, if the S(0-2) lines are major coolants, and if 
the line ratios are also consistent, a shock speed of order 
8 kms"-'^ would be sufficient to produce the observed line 
fluxes. Figure 4 of Draine et al. ( 1983) indicates that the 
H2 lines are in fact expected to dominate the cooling of 
Vg = 5 — 20 kms""'^ C-type shocks propagating into gas 
with Uh^o = 100 cm~^ and 7?o = 10 pG. 

Matching the line ratios is a separate challenge. The 
J = 3 and J = 4 levels could be excited in shocks. 
IDraine fc Katd (|1986f ) show results for shocks propagat- 
ing into gas with Uh = 50 cm"-^ and Bq = 7 pG. For 
Vg ~ 12 kms~^ and cos6' — 0.5, these models reproduce 
the observed S(2) emission and overestimate by a factor 
of about 2.6 the S(l) emission. The models assumed an 
ortho/para ratio of 3; using a smaller ratio (eg., 7 = 0.7) 
could decrease the S(l)/S(2) ratio to agree with the ob- 
served values. 

If such shocks are present in these clouds, one would ex- 
pect to see 21 cm emission with components correspond- 
ing to both the pre-shock and post-shock material. CO 
( J = 1 — 0) line emission should also show pre-shock and 
post-shock components. Such signatures have yet to be 
detected. 

Published shock models fail to predict the S(0) emis- 
sion we observe towards HLCs without overestimating 
the S(l) and/or S(2) lines. The very large column den- 
sity N{J = 2) — 1.5 X 10"'^^cm~2 appears to require a 
large filling factor of H2 gas with relatively cool tem- 
peratures. For thermal excitation we reproduce the S(0) 
emission if most of the gas has an excitation temperature 
of T = 85K (Fig. H black curve). 

One approach to increasing the J — 2 popula- 
tion relative to J = 3 and 4, while still reproduc- 
ing the overall column densities, has been to posit 
mult iple su perimposed, randomly oriented shock w aves. 
BothlGT^e l. Pineau des Forets. &: FedermanI (j2002l ) and 



Falgarone et~aIll2005D . 



Falgarone et al. (2005) applied the MHD shock model of 
Flower & Pineau des Forets ( 1998) to measurements of 
molecular lines in diffuse clouds and found that many 
of the observations could be explained by ensembles of 
shock waves. For comparison with our data, we dis- 
play on Fig. m (dash-dot curve) the sum of 8 shocks 
traveling with velocity Vg — S kms~^ into a medium of 
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nn = 50 cm~^ and B = 7 iiG ( models taken from Fig- 
ure 3 of iFalgarone at all |2005[) . iFalgarone et al.l found 
that, per magnitude of visual extinction, 3.2 shocks with 
Vs = 8 kms~^ can reproduce their S(l) and higher emis- 
sion from the diffuse medium. This is nearly the same 
"shock fraction" we need for the S(l) and S(2) data to- 
wards DCld 300.2-16.9 and LDN 1780 (our sample has 
« 2 mag, giving 4 shocks per magnitude). Since the 
shocks each have a column density iV^i, = 1.5 x 10^^ cm^^ 
of hot gas, the mass fraction of shocked gas in this model 
is about 0.06%. The model J — 3 and J = 4 column 
densities are reasonably well-matched to the data, but 
J = 2 is still insufhc i ently populated. Indeed, referring to 
the lDraine fc Katd (|1986l ) models, even a = 5 kms 
shock overproduces S(l) emission relative to S(0) by a 
factor of 8, so the J = 2 levels remain a problem to 
reproduce for shock models. 

Although shocks provide one mechanism for localized 
heating, it has also been proposed that intermittent dis- 
sipation of turbule nce can produ ce localized , transient 
pockets of hot gas ( Joulain et al.l [l998: Fale arone et al.l 
[20051: iGodardet al.ll2009[) . Weshow on Fig. H the ex- 
citat ion diagram for 1300 vortices, scaled from Figure 

3 of IFalgarone et al.l (|2005D (dot-dot-dash curve). The 
vortices have ion-neutral drift velocity vd = 3.5 kms~^, 
Uh — 30 cm~^, and _B = 10 fiG. The model has very sim- 
ilar H2 excitation to the multiple shock model described 
above, although it has a lower ortho/para ratio which 
gives a slightly better match to our S(l ) and S(2) data . 
We get the same "vortex fraction" as IFalgarone et al.l . 
~ 700 vortices per magnitude of visual extinction. The 
vortex column density of N^^r = 1.7 x 10 
mass fraction in the vortices of about 0.01% 

Since all cirrus is exposed to the ISRF, a more realis- 
tic model contains both UV and coUisional heating. We 
display on Fig. 2] the sum of the Gq = 1 FDR and the 
Vd = 3.5 kms""'^ vortex models (orange curve). The S(0) 
line remains underpredicted, by about a factor of 10. 

We have computed an additional model that repro- 
duces the H2 level populations quite well, even though it 
lacks a specific mechanism for heating the warm molec- 
ular gas (Draine & Ingalls, in preparation). This model 
is plotted in blue on Fig. S) It consists of a plane- 
parallel molecular cloud with an H i layer (FDR) on both 
front and back sides, superimposed with a zone of warm 
gas. The PDR portion of the slab obtains a uniform gas 
pressure P/k — 4000 cm"'^ K (consistent with pressure 
measurements in the ISM based on [C l] excitation — see 
iJenkins fc Tripp"b 011). and is embedded i n a radiation 
field equal to 0.6 times the iHabind ()1968[ ) energy den- 
sity at lOOOA. The warm zone is assumed to have been 
suddenly heated, so that the ortho/para ratio has not 
had time to reach equilibrium corresponding to T„^rm- 
CoUisional excitation of the H2 rotational levels is cal- 
culated using state-to-state rates. The data are best fit 
by r^arm = 391 K, an ortho/para column density ratio 
in both the warm and cool gas of 0.15, and a warm H2 
column density N„^^^(H.2) « 5.15 x lO^^cm"^, which is 
about 2% of the total H2 column density in the model. 
The bulk of the H2 emission in this model is not due to 
radiative pumping; more than 97% of the J = 2, 3, and 

4 excitation is produced in the warm component. 



cm ^ gives a 



4.3. II2 Excitation in the Supernova- Driven ISM 

In summary, the (somewhat speculative and incom- 
plete) scenario we envision is one in which CNM molec- 
ular clouds immersed in the interstellar radiation field 
of intensity Gq ^ 1 are subjected to a large scale in- 
jection of mechanical energy and dissipate that energy 
intermittently on smaller scales via shocks, or perhaps in 
localized transient turbulent dissipation regions (TDRs; 
Godard et al. 2009) such as vortices or velocity shears 
(or some combination thereof). The H2 S(0) - S(2) lines 
are presumed to be mainly excited by collisions in hot 
regions, heated by mechanical energy dissipation. The 
gas is initially cold (T < 100 K), and in the shielded in- 
terior of a cloud the ortho/para ratio may be very low. 
When dissipational heating occurs, local temperatures 
temporarily reach values near 1000 K, and ortho/para 
rises. If the source of large scale mechanical energy is 
not refreshed somehow, the fraction of warm gas will de- 
crease until the line radiation is no longer detectable. 

Since both clouds in this study are thought to have 
been affected by supernova explosions from t he Sco-Cen 
OB association about (2 — 6) x 10^ yr ago (jToth et al.l 
[l99l iNeh^e et al.ll2Q08ah . we surmise that portions of 
the medium are still dissipating the energy deposited 
from those 100 pc scale blast waves. For this scenario 
to hold, the time to dissipate the gas turbulent kinetic 
energy must be larger than the time since blast wave im- 
pact. Using the total II2 S(0 — 2) intensity from the DCld 
300.2-16.9 positions, /(0-2) w 5 x 10"^ Wm^^ sr^i, we 
can calculate the H2 cooling rate per H atom: 



Ah- 



47r/(0- 2) 

3.1 X 10"^^W(H atom)-^ 

/(0-2) \ /2 X 10^1 cm-2 



5 X 10-9Wm-2sr-i 
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(11) 



The turbulent kinetic energy per proton is 



1.4 X — m„ (T, 



3.2 X 10"^*^ J (H atom)-i f , , 

\3 kms ^ 



(12) 



where mp is the proton mass, Uv is the gaussian standard 
deviation of the observed velocity distribution, and the 
factor 1.4 accou nts for helium. Usin g Parkes Galactic 
All Sky Survey (jKalberla et al.ll20ldf ) H i observations 
towards DCLd 300.2-16.9, we measure k, Z kms~^. 
Dividing the turbulent energy by the cooling rate gives 
the dissipation time: 



^x^«^-(3k^ 
2 X 1021 cm-2 



2 /5 X 10-9Wm-2sr-i 



/(0-2) 



(13) 



Thus, the time over which the observed II2 emission can 
dissipate the available turbulent kinetic energy is of the 
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same order of magnitude as the time since supernova 
impact. 

This dissipation time might be an upper hmit, since 
the H2 hues are not the only coolant of the warm gas. 
Above about 200 K the 63/im fine-structure line of neu- 
tral oxygen [O i] can account for a significant fraction of 
the cooling; and below 200 K the 158/i m line of ionized 
carbo n [C li] can dominate the cooling (iFalgarone et al.l 
[2OOI . Our non-LTE+PDR model predicts 63/xm [O i] 
and 158/xm [C 11] emission from the T — 391 K warm zone 
of « 1 X 10~^ Wm^^ sr^^ each, which increases the total 
cooling by about 40% over that due to H2 emission alone, 
leading to a corresponding decrease in the estimated dis- 
sipation time. (The model assumes an oxygen abundance 
of 4 X 10~^ and a carbon abundance of 1 x 10"''.) For 
DCld 300.2-16.9, the reduced dissipation time is about 
the same as the time si nce supernova impact deduced 
by iNehme et all (|2008aD (2 - 3 x 10^ yr, based on the 
time it took the observed 10 kms~^ H i streaming mo- 
tions to create the cloud's 2.4 pc tail). The detection 
of [O i] 63/im emission at the ~ 10~^Wm~^sr~^ level 
would help to confirm our assertion that the H2 emis- 
sion we observe originates in anomalously warm regions, 
since the model predicts about 6 times as much [O i] 
flux in the warm zone as in the PDR. Nearly all of the 
HLC positions in o ur study show [C 11] 158/im emission 
(jlngalls et al.l 120021 ) at a level consistent with the PDR 
part of the model. Since only 5% of [C 11] emission is ex- 
pected to originate in the 391 K zone, this cannot likewise 
be used to corroborate the existence of warm material. 

For LDN 1780, we have inferred a time of 6 x 10^ yr 
since supernova impact, which is somewhat larger than 
the estimated dissipation time. We lack S(0) or S(l) 
measurements towards this cloud, however, so it is not 
clear that <dis estimated from the DCld 300.2-16.9 total 
S(0-2) intensity is appropriate for LDN 1780. 

The timescale to expend the available kinetic energy 
can also increase by a factor of a few or more, depending 
on the specific mechanism of dissipation. If the H2 emis- 
sion is primarily the result of large-scale shocks, then 
dissipation will begin soon after energy injection. If 
the emission is caused by turbulent dissipation at small 
scales, however, then any large-scale mechanical energy 
input will take additional time to flow from the input 
scale to the dissipation scale. For a Kolmogorov-like tur- 
bulent spectrum the dissipationless cascade will be dom- 
inated by the large-scale turnover time, L/<7„. Given a 
velocity dispersion of 3 kms~^ and an input size scale of 
L ^ 1 pc, this turnover time is « 3 x 10*^ yr, of order the 
dissipation time itself. 

In general, a given spot in the Galactic disk experi- 
ences a supernova explosion within 100 pc once every 
~ 10^ yr. Since any energy deposited would dissipate 
in idis ~ 3 X 10^ yr, our scenario provides a natural ex- 
planation for the statistical dearth of rotationally-excited 
H2 in other translucent regions in our 100 pc neighbor- 
hood (we only detected 12.3/im emission towards 2 out 
of 6 HLCs observed). The case of LDN 183, consid- 
ered part of the sarne molecular cornplex a s LDN 1780 
(jLaureiis et al.l [19951 iLallement et al.ll2003[ ). lends plau- 
sibility to this picture. LDN 183, which lacks detectable 
S(2) emission, is outside the arc of the H i loop that 
passed through LDN 1780 ()T6th et al.lll995l) ~ 6 x 10^ yr 



ago, and so perhaps has not yet been subjected to the 
same large scale energy injection that excited the S(2) 
emission we see towards LDN 1780. 

4.4. A New Source 0/II2 Photons from Galaxies? 

DCld 300.2-16.9 and LDN 1780 could be templates 
for an as-yet undefined "translucent" source of H2 pure- 
rotational emission from the disks o f all non-active galax- 
ies. For a large sample of galaxies, iRoussel et al.l |2007D 
found that the ratio between the total H2 flux in the 
S(0)-S(2) transitions and the PAH emission at 7.9/im, 
F(S0-S2)/F7 PAH, was remarkably constant. The ratio 
was constant regardless of whether the IR emission was 
dominated by cirrus or high-intensity PDRs (their Fig- 
ure 9), which they took to signify that the emission from 
both PAHs and H2 is caused chiefly by FUV excitation in 
PDRs. But when we compare our Sp itzer F(S2)/F7_9 pah 
measurements to the IRoussel et al.l data (Figure [5]), we 
see that cirrus positions are paradoxically brighter in S(2) 
emission (relative to PAH emission), by a factor of 4 
(on average). Thus translucent material appears to be 
at least as effective at producing rotationally excited H2 
molecules (per PAH-exciting photon) as PDRs surround- 
ing H II regions. 

Including the S(0) and S(l) emission for DCld 300.2- 
16.9 and estimating the TIR emission yields S(0-2)/TIR 
values that are ~ 10 times higher on average than the 
non-active galaxy ensemble. This implies that the over- 
all H2 excitation in translucent cirrus is higher per unit 
of total energy absorbed by dust from the ambient ra- 
diation flel d. Qualitatively, this is not a controversial 
statement; iBlack and van DishoecU (|1987[ ) showed that 
regions with Go/n <C 1 are more efficient at radiatively 
exciting H2 than regions with Go/n > 1, because when 
Go /n is high the resulting low H2 abundance means that 
the IIIO-912A UV photons that might pump H2 are in- 
stead absorbed by dust grains. But as we have shown, 
the S(0) to S(2) emission we observe is far greater than 
that expected from PDR models, requiring us to invoke 
collisional excitation in anomalously warm regions. 

The HLCs in which we have detected S(2) emis- 
sion have a greater concentration of rotationally warm 
gas than the integrated disks of galaxies — including the 
Milky Way. The integrated S(2) emission from a col- 
umn of material equivalent to ~ 30 mag of visu al extinc- 
tion in the Milky Way (F algarone et al.l 120051 ) is indis- 
tinguishable on Figure [S] from the SINGS star-forming 
galaxies. One-half of the mass in the MW column i s 
calculated to be in diffuse gas ([Falgarone et al.l [20051) . 
so it is reasonable to expect that other galaxies likewise 
have a significant mass of diffuse gas that emits S(2) 
radiation. Indeed, there is evidence that cold gas with 
detectable S(2) emission, if not common, may neverthe- 
less be present throughout the disks of normal galax- 
ies. As mentioned, the data are consistent with other 
measurements of nearby molecular lin es of sight: anal- 
ysis o f both the emission from Taurus (!Gold smith et al.l 
|2010( ) and local absorption towards stars me asured with 
FUSE fR achford. et al. 11200 It IgTV et al.ll200l shows sim- 
ilar excitation conditions to those reported here (§3.2). 
In addition, recent measurements of H2 in nearby galax- 
ies (|Laine et al.ll20Tot IStacev et al.ll2010l) yield S(0)/S(1) 
rotational temperatures that are similar to those we niea- 
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sure (~ 100 — 150 K). These temperatures do not vary 
across the disks of the galaxies, despite the fact that star 
formation indicators decrease radially, suggesting a com- 
ponent of excitation that is unrelated to PDRs surround- 
ing O and B stars. 

Assume that the 34-position high-latitude cloud sam- 
ple is representative of all translucent gas in the Galaxy. 
If so, then the average F(S2)/F7 9_pah ratio for all po- 
sitions, including nondetections, estimates the overall 
Galactic ratio for such material. A gray dashed ellipse 
on Figure [S] represents this average, and is centered on 
practically the same value of F(S2)/F7 9 pah as the SINGS 
average. It is therefore plausible that translucent clouds 
integrated over the disks of galaxies represent a signifi- 
cant source of the II2 pure-rotational emission. 

4.5. Future Work to Characterize Warm H2 in Cirrus 

Clouds 

We have proposed a rather speculative scenario to ex- 
plain the S(0) through S(2) emission we detect towards 
two translucent cirrus clouds. We used models of radia- 
tive and mechanical heating to estimate the conditions 
that could reproduce the observations. Knowledge of the 
superbubble environment of the clouds lends credence to 
the idea that extra mechanical energy is available to raise 
enough H2 to the J = 4 state. The processes responsible 
for the heating remain uncertain. T urbulent dissipation 
regions (TDRs; IGodard et ai1[200l may be able to pro- 
duce strong localized heating, but more complete theo- 
retical models of how this may proceed in magnetized 
interstellar gas remain to be developed. If shock waves 
are responsible for the heating, it should be possible to 
identify other evidence for the presence of a shock, such 
as separate radial velocities in the preshock and post- 
shock regions. 

New observations are also required to confirm our in- 
terpretation of the data, especially our assertion that H2 
rotational emission from cirrus is widespread in normal 
galaxies, and that cirrus gas is a major source of the inte- 
grated emission. H2 absorption line studies with FUSE 
have already established the ubiquity of H2 in J = 2, 
3, 4, and 5 at levels in excess of what can be explained 
by UV pumping ( Wakker ^ 2006; Gillmon et al. 2006) — 
see Figure 31.4 in iDraind (j201l[) . The Cosmic Origins 
Spectrograph on Hubble Space Telescope can make ad- 
ditional H2 absorption line measurements in diffuse and 
translucent clouds when suitable background UV sources 
are available. Later in this decade, it will be possi- 
ble to measure the IR line emission from rotationally- 
excited H2 using both the James Webb Space Telescope 
(JWST) and the proposed Space Infrared Telescope for 
Cosmology and Astrophysics (SPICA). Other tracers of 
warm diffuse molecular gas, such as CH+ a nd SH+ (e.g., 
iMenten et aD[20Tlt iFaJgarone et al.ll2010al|0 ). and [O i], 
can also be used to study the heating processes in diffuse 
and translucent material. 

5. CONCLUSIONS 

We have detected molecular hydrogen pure-rotational 
(u = — 0) emission up to the S(2) transition towards 
6 positions in two nearby high Galactic latitude translu- 
cent cirrus clouds (HLCs). Excitation analysis for the 
three available lines show that a small amount (< 2%) 



of the gas is rotationally warm (T^ > 300 K). Pho- 
todissociation region models with Go = 0.6 — 10 are un- 
able to account for the observed emission, whereas mod- 
els of shocks or turbulent dissipation regions can repro- 
duce the S(l) and S(2) measurements. Both clouds have 
likely been impacted by shock waves from the Scorpius- 
Centaurus OB association about 2 — 6 x 10^ yr ago, which 
may have injected sufficient mechanical energy to excite 
the H2 rotational levels. Four other nearby HLCs, which 
are not so situated, did not show detectable S(2) lines. 
This picture is still speculative, however, and requires 
robust modeling to confirm it. 

Prior to the Spitzer Space Telescope, published obser- 
vations of molecular hydrogen in the diffuse interstellar 
medium were confined to UV absorption measurements 
(for example, se c Rachford, ct al. 2001, 2002; Grv ct aO 
12002; Gillmon ct al. 2006; Wakkcr 2006; Rachford. et all 
2009; Jensen et aL .2010) or a measurement of the emis- 
sion along an extensive ( Ay ^ 30 mag) line o f sight 
in the Galactic midplane (jFalgarone et al.l [20051 ). Dur- 
ing its cryogenic lifetime, Spitzer has enabled a new 
way of looking at the neutral ISM, confirming the no- 
tion that a significant amount of warm molecular gas 
exists within the cold medium, both w ithin our own 
Galaxy (this paper: iGoldsmith et al.ll2010D. a nd in exter- 
nal galaxies (jLaine et al.l 120101: iStacev et al.ll2010[ ). We 
propose that rotationally warm gas inside cold clouds 
is a feature of the disks of all normal galaxies. While 
it may not be common (less than 20% of our sample 
had detectable H2 emission, and at most 2% of the 
column density towards the detected positions is warm 
gas), it appears to be widespread. The enhanced H2 
flux relative to dust emission in HLCs compared with 
the inner disks of the SINGS galaxies, together with the 
known ubiquity of ^ 1 material in the Galactic plane 
(|Polk. Knapp. Stark, fc WilsonI flQSSh . requires that the 
translucent cirrus be taken into account when interpret- 
ing the H2 rotational lines in the disks of all non-active 
galaxies. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. Support for this work 
was provided by NASA through an award issued by 
JPL/Caltech. We thank the anonymous referee for com- 
ments which improved this paper. Facilities: Spitzer () 
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